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Abstract 

We review some recent developments which make use of the con- 
cept of 'superstatistics', an effective description for nonequilibrium 
systems with a varying intensive parameter such as the inverse tem- 
perature. We describe how the asymptotic decay of stationary prob- 
ability densities can be determined using a variational principle, and 
present some new results on the typical behaviour of correlation func- 
tions in dynamical superstatistical models. We briefly describe some 
recent applications of the superstatistics concept in hydrodynamics, 
astrophysics, and finance. 



1 Introduction 



Complex nonequilibrium systems often exhibit dynamical behaviour that is 
characterized by spatio-temporal fluctuations of an intensive parameter f3. 
This intensive parameter may be the inverse temperature, or an effective 
friction constant, or the amplitude of Gaussian white noise, or the energy 
dissipation in turbulent flows, or simply a local variance parameter extracted 
from a signal. A nonhomogeneous spatially extended system with fluctua- 
tions in (3 can be thought of as consisting of a partition of spatial cells with 
a given f3 in each cell. If there is local equilibrium in each cell (so that sta- 
tistical mechanics can be applied locally), and if the fluctuations of (3 evolve 
on a sufficiently large time scale, then in the long-term run the entire system 
is described by a superposition of different Boltzmann factors with different 
(3, or in short, a 'superstatistics' [T|. The superstatistics approach has been 
the subject of various recent papers I^ElSlfeliniEZlIHlO HOI HH H21 UHl - Su- 
perstatistical techniques can be successfully applied to a variety of physical 
problems, such as Lagrangian fHl ITo] and Eulerian turbulence [THl HZ], defect 
turbulence [TH], cosmic ray statistics [TH], plasmas PU]. statistics of wind veloc- 
ity differences |21[ |2*2~] and mathematical finance j2Hl I2U I2H]- Experimentally 
measured non-Gaussian stationary distributions with 'fat tails' can often be 
successfully described by simple models that exhibit a superstatistical spatio- 
temporal dynamics. 

If the intensive parameter in the various cells is distributed according to a 
particular probability distribution, the ^-distribution, then the correspond- 
ing superstatistics, obtained by integrating over all (3, is given by Tsallis 
statistics 123 HE1 For other distributions of the intensive parame- 
ter (3, one ends up with more general superstatistics, which contain Tsallis 
statistics as a special case. Generalized entropies (analogues of the Tsallis 
entropies) can also be defined for general superstatistics [HJ E2J EH EH] and are 
indeed a useful tool. The ultimate goal is to proceed from simple models to 
more general versions of statistical mechanics, which are applicable to wide 
classes of complex nonequilibrium systems, thus further generalizing Tsallis' 
original ideas [26 . 

This paper is organized as follows: 

First, we briefly review the superstatistics concept. We then show how 
one can deduce the asymptotic decay rate of the stationary probability den- 
sities of general superstatistics from a variational principle [3Tj. In section 4 
we consider dynamical realization of superstatistics and investigate the typ- 
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ical behaviour of correlation functions. Many different types of decays of 
correlations (e.g. power law, stretched exponentials) are possible. Section 5 
summarizes some recent applications of the superstatistics concept in hydro- 
dynamics, astrophysics, and finance. 

2 What is superstatistics? 

The superstatistics approach is applicable to a large variety of driven nonequi- 
librium systems with spatio-temporal fluctuations of an intensive parameter 
/3, for example, the inverse temperature. Locally, i.e. in spatial regions (cells) 
where (3 is approximately constant, the system is described by ordinary statis- 
tical mechanics, i.e. ordinary Boltzmann factors e~@ E , where E is an effective 
energy in each cell. In the long-term run, the system is described by a spatio- 
temporal average of various Boltzmann factors with different (3. One may 
define an effective Boltzmann factor B(E) as 

POO 

B(E)= / f((3)e^ E d(3=(e^ E } } (1) 
Jo 

where /(/?) is the probability distribution of (3 in the various cells. For 
so-called type-A superstatistics lj, one normalizes this effective Boltzmann 
factor and obtains the stationary long-term probability distribution 

P(E) = \B{E\ (2) 

where 

Z= B(E)dE. (3) 
Jo 

For type-B superstatistics, the /3-dependent normalization constant of each 
local Boltzmann factor is included into the averaging process. In this case 
the invariant long-term distribution is given by 

p(£) =i m m dP ' (4) 

where Z((3) is the normalization constant of e~ l3E for a given f3. Eq. (J3J) 
is just a simple consequence of calculating marginal distributions. Type- 
B superstatistics can easily be mapped into type-A superstatistics by re- 
defining /(/?). 
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A superstatistics can be dynamically realized by considering Langevin 
equations whose parameters fluctuate on a relatively large time scale (see 
[3*2] for details). For example, for turbulence applications one may consider a 
superstatistical extension of the Sawford model of Lagrangian turbulence [HJ 
This model consists of suitable stochastic differential equations for 
the position, velocity and acceleration of a Lagrangian test particle in the tur- 
bulent flow, and the parameters of this model then become random variables 
as well. Experimental data are well reproduced by these types of models. 

Often, a superstatistics just consists of a superposition of Gaussian dis- 
tributions with varying variance. The parameter (3 can then be estimated 
from an experimentally measured signal u(t) as 



where (...)t denotes an average over a finite time interval T of the signal, 
corresponding to the 'cell size' of the superstatistics. It is then easy to make 
histograms of (3 and thus empirically determine f(/3). 

3 Asymptotic behaviour for large energies 

Superstatistical invariant densities, as given by eq.Q or (J1J), typically exhibit 
'fat tails' for large E, but what is the precise functional form of this large 
energy behaviour? The answer depends on the distribution f{0) and can 
be obtained from a variational principle. Details are described in |3T], here 
we just summarize some results. For large E we may use the saddle point 
approximation and write 
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(5) 






(6) 



where 



0s = sup{-/3£ + In /(/?)}. 



(7) 
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The expression 

sn P {~PE + In f((3)} (8) 

P 

corresponds to a Legendre transform of ln/(/3). The result of this transform 
is a function of E which can be thought of as representing a kind of entropy 
function if we consider the function In f{(3) to represent a free energy function. 
This entropy function, however, is different from other entropy functions used 
e.g. in nonextensive statistical mechanics. It describes properties related to 
the fluctuations of inverse temperature. 

In the case where f(/3) is smooth and has only a single maximum we can 
obtain the supremum by differentiating, i.e. 

sup{-/3£ + In /(/?)} = -faE + ln/(/fe) (9) 

P 

where (3e satisfies the differential equation 

= -E + (\nf(f3)Y = -E + ^. (10) 

By taking into account the next-order contributions around the maximum, 
eq. © can be improved to 

\/-(ln /(&))" 

Let us consider a few examples. Consider an f(/3) of the power- law form 
f(j3) ~ /i 7 , 7 > for small (3. An example is a \ 2 distribution of n degrees 
of freedomjnaiSl!, 

'<*>=4te) (12) 

(A) > 0, n > 1) which behaves for (3 — > as 

/(/3) ~ /3"/ 2 - 1 , (13) 

7 = | - 1- (14) 
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Other examples exhibiting this power-law form are F-distributionsPU l2*Uj . 
With the above formalism one obtains from eq. (jl(Jj) 

= | (15) 

and 

B{E) ~ E" 7 - 1 . (16) 

These types of f{(3) form the basis for power-law generalized Boltzmann 
factors (g-exponentials) B(E), with the relation H7\ l2l ^ 123] 

7+1 = ^- (17) 

Another example would be an f(/3) which for small j3 behaves as f{(3) ~ 
e -c//3^ c > g j n ^ g cage Qne OD t a i n s 

i>. = yjl (18) 

and 

The above example can be generalized to stretched exponentials: For f(/3) 
of the form f{(3) ~ e~ c/3 one obtains after a short calculation 



and 

B^-E^^e^' 1 ', (21) 

where a is some factor depending on 5 and c. 

Of course which type of f(/3) is relevant depends on the physical system 
under consideration. For many problems in hydrodynamic turbulence, log- 
normal superstatistics seems to be working as a rather good approximation. 
In this case /(/?) is given by 

hr\ l / (log/3-'") 2 ! ,,,, 

m = 7srr p \ — w— (22) 



where s and m are parameters [TJ HH HH1 HEl UH EE] • 
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4 Superstatistical correlation functions 



To obtain statements on correlation functions, one has to postulate a concrete 
dynamics that generates the superstatistical distributions. The simplest dy- 
namical model of this kind is a Langevin equation with parameters that vary 
on a long time scale, as introduced in [32]- 

Let us consider a Brownian particle of mass m and a Langevin equation 
of the form 

v = - 1V + aL(t), (23) 

where v denotes the velocity of the particle, and L(t) is normalized Gaussian 
white noise with the following expectations: 

(L(t)) = (24) 
(L(t)L(t')) = 5(t-t'). (25) 

We assume that the parameters a and 7 are constant for a sufficiently long 
time scale T, and then change to new values, either by an explicit time 
dependence, or by a change of the environment through which the Brownian 
particle moves. Formal identification with local equilibrium states in the cells 
(ordinary statistical mechanics at temperature /5 _1 ) yields during the time 
scale T the relation [3^] 

(v 2 ) = — = -^- (26) 
x 1 2 7 (3m K ' 

or 

2 7 

P = -~r (27) 
m a z 

Again, we emphasize that after the time scale T, 7 and a will take on new 
values. During the time interval T, the probability density P(v,t) obeys the 
Fokker-Planck equation 



dP d(vP) 1 2 d 2 P 
dt ^ dv 2° dv 



7^ + -y"^ (28) 



with the local stationary solution 



PH/?) = y^eKpi-i/Wj. (29) 

In the adiabatic approximation, valid for large T, one asumes that the local 
equilibrium state is reached very fast so that relaxation processes can be 
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neglected. Within a cell in local equilibrium the correlation function is given 
by [M| 

C{t-t'\(3) = (v(t)v{t')) = ^-e-^- t \ (30) 
Clearly, for t = t' and setting m = 1 we have 

p= (k' (31) 

in agreement with eq. . 

It is now interesting to see that the long-term invariant distribution P{v), 
given by 

POO 

P(v) = / f(ft)P(v\ft)dft (32) 

JO 

depends only on the probability distribution of ft = —\ and not on that 
of the single quantities 7 and a 2 . This means, one can obtain the same 
stationary distribution from different dynamical models based on a Langevin 
equation with fluctuating parameters. Either 7 may fluctuate, and a 2 is 
constant, or the other way round. On the other hand, the superstatistical 
correlation function 

poo -1 poo 

C(t -f)= f(P)C(t - t'\(3)d(3 = - / f^p-'e-^dp (33) 
Jo m Jo 

can distinguish between these two cases. The study of correlation functions 
thus yields more information for any superstatistical model. 

Let illustrate this with a simple example. Assume that a fluctuates and 
7 is constant such that f3 = — is x 2 -distributed. Since 7 is constant, we 
can get the exponential e~ 7 ' t_t ' out of the integral in eq. (J3*3*j) . meaning that 
the superstatistical correlation function still decays in an exponential way: 

C{t-f) -e- 71 *-*' 1 . (34) 

On the other hand, if a is constant and 7 fluctuates and ft is still x 2 - 
distributed with degree n, we get a completely different answer. In this case, 
in the adiabatic approximation, the integration over ft yields a power-law 
decay of C{t-t'): 

C{t-t')~\t-t'\-\ (35) 
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where 

n 

V = g " 1 ( 36 ) 

Note that this decay rate is different from the asymptotic power law decay 
rate of the invariant density P(v), which, using (|2TIJ) and (JH2j), is given by 
P(v) ~ v-V^-V, with 

1 n 1 . 

+ o- (37) 



g- 1 2 2 

In general, we may generate many different types of correlation functions 
for general choices of f{(3). By letting both a and 7 fluctuate we can also 
construct intermediate cases between the exponential decay (JMj) and the 
power law decay ()35|) . so that strictly speaking we only have the inequality 

V > g - 1. (38) 

depending on the type of parameter fluctuations considered. 
One may also proceed to the position 



x {t) = I v {t')dt' (39) 

of the test particle. One has 







( x \t)) = [ [\v(t')v(t"))dt'dt". (40) 
Jo Jo 

Thus asymptotic power-law velocity correlations with an exponent 77 < 1 are 
expected to imply asymptotically anomalous diffusion of the form 

(At)} - f (41) 

with 

a = 2-r]. (42) 

This relation simply results from the two time integrations. 

It is interesting to compare our model with other dynamical models gener- 
ating Tsallis statistics. Plastino and PlastinojnZl and Tsallis and Bukmann 38J 
study a generalized Fokker-Planck equation of the form 

Em = -£ (Hx)Plx , t))+D <jL P(x , t r (43) 
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with a linear force F(x) = k% — k 2 x and v ^ 1. Basically this model means 
that the diffusion constant becomes dependent on the probability density. 
The probability densities generated by eq. (JlHj) are g-exponentials with the 
exponent 

q = 2-u. (44) 

The model generates anomalous diffusion with a = 2/(3 — q). Assuming the 
validity of a = 2 — f), i.e. the generation of anomalous diffusion by slowly 
decaying velocity correlations with exponent fj, one obtains 

On the other hand, for the x 2 -superstatistical Langevin model one obtains 
by combining eq. (|3*Bj) and (}3Tj) the different relation 

5 - 3g 

71 = 2q^2- (46) 

Interesting enough, there is a distinguished g-value where both models yield 
the same answer: 

q = 1.453 rj = rj = 0.707 (47) 

These values of q and rj correspond to realistic, experimentally observed 
numbers, for example in defect turbulence |18j. 

So far we mainly studied correlation functions with power law behaviour. 
But in fact one can construct superstatistical Langevin models that exhibit 
more complicated types of asymptotic behaviour of the correlation functions. 
To see this we notice that the asymptotic analysis of section 3 applies to 
correlation functions as well, by formally defining 

E := ^a 2 m\t-t'\ (48) 

/(/?) = (49) 



and writing 



oo 



C(t~t')= / f(P)e- pE dp. (50) 



To obtain statements on the symptotic decay rate of the superstatistical 
correlation function, we may just use the same techniques described in section 
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3 with the replacement E — > E and f f. In this way one can construct 
models that have, for example, stretched exponential asymptotic decays of 
correlations etc. (see also 39 ). Asymptotic means here that \t — 1'\ is large 
as compared to the local equilibrium relaxation time scale, but still smaller 
than the superstatistical time scale T, such that the adiabatic approximation 
is valid. 

5 Some Applications 

We end this paper by briefly mentioning some recent applications of the su- 
perstatistics concept. Rizzo and Rapisarda[2Tl 122] study experimental data 
of wind velocities at Florence airport and find that x 2 -superstatistics does 
a good job. Jung and Swinnev[33] study velocity differences in a turbu- 
lent Taylor- Couette flow, which is well described by lognormal superstatis- 
tics. They also find a simple scaling relation between the superstatistical 
parameter (3 and the fluctuating energy dissipation e. Paczuski et al. 40 
study data of solar flares on various time scales and embedd this into a 
superstatistical model based on x 2 -superstatistics = Tsallis statistics. Hu- 
man behaviour when sending off print jobs might also stand in connection 
to such a superstatistics [H]. Bodenschatz et al.[I2] have detailed experi- 
mental data on the acceleration of a single test particle in a turbulent flow, 
which is well described by lognormal superstatistics, with a Reynolds num- 
ber dependence as derived in a superstatistical Lagrangian turbulence model 
studied by Reynolds [To']. The statistics of cosmic rays is well described by 
X 2 -superstatistics, with n = 3 due to the three spatial dimensions[19j. In 
mathematical finance superstatistical techniques are well known and come 
under the heading 'volatility fluctuations', see e.g.[23j for a nice introduction 
and [2U |25] for some more recent work. Possible applications also include 
granular media, which could be described by different types of superstatistics, 
depending on the boundary conditions [43]. The observed generalized Tsal- 
lis statistics of solar wind speed fluctuations^!] is a further candidate for 
a superstatistical model. Chavanis[3Ti] points out analogies between super- 
statistics and the theory of violent relaxation for collisionless stellar systems. 
Most superstatistical models assume that the superstatistical time scale T 
is very large, so that a quasi-adiabatic approach is valid, but Luczka and 
Zaborek^3] have also studied a simple model of dichotomous fluctuations of 
j3 where everything can be calculated for finite time scales T as well. 
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